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We investigate the fate of bottom quarks produced in heavy-ion collisions at RHIC. 
Examining both the direct capture of a charmed quark, and multi-step processes, 
where the B c meson is formed in a sequence of quark capture and exchange re- 
actions, we find: a) that a sufficiently high number of B c 's will be produced to 
generate a detectable tri-lepton signal, and b) that the production rate of B c 's is 
highly sensitive to the properties of the deconfined source. A flavor-independent 
potential model, which includes color screening effects, is used to study the prop- 
agation of a B c in a quark-gluon fireball and to compare this behavior with that 
of J/ip mesons. 

1 Introduction 

We argue here that RHIC energies should create physical conditions conducive 
to the use of B c mesons as a new observable of deconfmement and quark-gluon 
plasma (QGP) formation, in spite, and also .because of the fact that elementary 
processes lead to the small branching ratio El, 

Re iem = Bc+_B* C _ 10 _4_ 10 -b _ (1) 
bb 

The fast electromagnetic decay of the vector meson B* allows us to ignore here 
the distinction between the scalar and vector mesons. The CDF collaboration 
recently announced the discovery of the B c systemBtl. Their announcement was 
based on careful analysis of 20 events with characteristic three-charged particle 
tracks expected from the decay B^r — ► Jji\) + P 1 + X. The CDF collaboration 
determined values for both the mass and the lifetime of the ground state, 

M(B C ) = 6400 ± 390 ± 130 MeV, t{B c ) = 0.46±g;J| ± 0.03 ps. (2) 

Allowing that either the c quark or the b quark can decay while the other 
is simply a spectator, or that the c quark and the b can annihilate into an 
intermediate vector boson, the weak decay width of the ground state of the B c 
is the sum of three terms, namely, 

T(B C -+X) = T(b -+X) + T(c -> X) + r(annih). (3) 
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Thus our potential model calculation of section || yields t(B c ) = 0.36 ± 0.05 ps, 
which agrees with the CDF result. Our calculation of the mass of B c described 
below in section || is also in good agreement with the CDF measurement. 

This lifetime allows B c mesons with sufficient transverse velocity to escape 
the beam centroid and thus produce a secondary tri-lepton vertex, which serves 
as a distinctive feature for identifying the decay. It is noteworthy that the life- 
times of all B c states below threshold are much longer than their counterparts 
in charmonium or the upsilon system because flavor conservation prevents de- 
cays through gluon annihilation diagrams. 

In a QGP environment at RHIC the formation of this 'stable' exotic B c 
meson is facilitated by the ready availability of charmed quarks. We explore 
here the production enhancement of B c with a view at both a study of its 
properties and as diagnostic tool of the QGP state. 

2 B c Production and Survival in a Nuclear Fireball 

Although the energy available to the individual nucleons in a typical RHIC 
experiment is not as high as that at Fermilab, the formation of QGP presents 
an environment that should enhance the production of bound B c states, since 
both b and c quarks should be able to propagate freely. At typical RHIC 
energies parton fusions, 

9 + 9 -> Q + Q, q + q^Q + Q, (4) 

should produce a supply of propagating charmed and occasional b and hjquarks 
with energies less than 10 GeV. The recent calculation of Mustafa et alu shows 
that gluon bremsstrahlung provides a very efficient mechanism for energy loss 
and that most of the heavy quarks with energies in this range should stop in 
a typical distance of 2 fm. An adequate flux of charmed quarks would lead to 
the production of B c states by an inelastic scattering accompanied by gluon 
bremsstrahlung. 

Since the strange quark flux is probably about 100 times larger than the 
charmed quark flux, B s states should be formed much more frequently than 
B c states. Because of the tight binding of the B c states, quark exchange 
B s +c <-> B c + s is an exothermic reaction, and we expect this two-step process 
to be also an important mechanism for production of -B c 's. 

Using a total bb production cross section of 1.7 microbarns and an overlap 
function appropriate for a head-on Au-Au collision at 100 A GeV, we get a 
bb production rate of 0.05 pairs per central collision. The abundance of cc- 
pairs has been estimated to be 200 times bigger, and our population evolution 
calculations yield a much greater fraction Rb c than has been obtained for the 
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one step reactions, Eq. (jl]). A further enhancement of the relative i3 c -yield 
arises from multi-step processes described above, operating in particular in 
the final stages of the QGP evolution. Our calculations show in presence of 
plasma a major change in the value Rb c , Eq. (0), and thus a major increase in 
the absolute yield of the B c mesons. 
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Figure 1: Quarkonium dissociation rates as functions of temperature. 



We cannot describe here all the intricacies of the population evolution 
equations that establish the final production rates of the B c . We mentioned 
already that there are several channels of production. Beyond that the issue 
is to assess the prospects for the survival of the B c state in the hot parton 
environment. Using Kharzeev's approach!], we have calculated the quarkonium 
dissociation rates for J/ip, B c and T(IS'), as shown in Fig.[|. The curves there 
indicate that B c should travel 5 to 10 times further than J/tp at T = 300 MeV. 
It is worth noting that the probability of dissociation of the B c at T = 400 MeV 
is about the same as that of J/ip at T — 300 MeV. 

3 Potential Model Approaches to the B c System 

Another way to understand the unusual stability of the B c mesons in the 
plasma phase is to consider its bound state structure. The rich spectra of 
bound states of charmonium and the upsilon system below the flavor thresh- 
olds □ provide enough information to determine the parameters of a nonrela- 
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tivistic potential model for heavy quarks. The assumption of flavor indepen- 
dence extends this model to the B c system. We describe here one such recent 
calculation^! incorporating a model of a running coupling constant effects in the 
central potential and the full radiative one-loop expressions supplemented by 
the Gromes consistency condition to incorporate non-perturbative effects in the 
spin-dependent potentials: these calculations provided an excellent fit of the 
upsilon levels (avg. dev. = 4.3 MeV), a good fit of the charmonium levels (avg. 
dev. = f 9. 9 MeV) and a good account of the leptonic widths below thresh- 
old. The model predicts a ground state energy of M{B C ) — 6286tg 5 MeV and 
that this state is 820 MeV below the threshold for heavy-light flavored meson 
production. Notably, the ground state of B c is expected to be about 150 MeV 
more tightly bound than the ground state of charmonium, and one thus would 
expect it to be more likely to survive QGP environments with temperatures 
around 300 MeV. 

To extend the potential model concept to QGP, one needs a means to 
consider the effects of color screening. Thus, we follow the the approach of 
Karsch, Mehr and Satz 0, who introduce color screening effects with an ex- 
ponential damping factor. Their analysis begins with a Coulomb plus linear 
potential, which describes the T = limit. Such a potential is known to give 
an adequate account of heavy-quark spin-averaged energies Their inverse 
Debye screening length parameter /i(T) describes the temperature dependence 
of color screening, ft also allows for deconfinement by parameterizing the melt- 
ing of the string tension. The Karsch- Mchr-Satz form for the central potential 
is 

v KMS (p(T),R) = A {1 ~ e " T) -«— . (5) 

From the large distance limit, the energy required to dissociate the system is 

Ecti s ((JL) = m b + m c H E bc (fj,) , (6) 

where the binding energy £Jj, c (/z) is obtained solving the Schrodinger equation 
with Vkms(m(7 1 ) j R), Eq. (||). Deconfinement is achieved when the dissociation 
energy vanishes, thus defining the critical value of /i. Our results for the 
dissociation energies of the three lowest B c states are shown in the top portion 
of Fig. § There one sees that fi crit (lS) = 840 MeV, almost 200 MeV higher 
than the corresponding quantity for J/ip. The [i dependence of the sizes of 
these 3 lowest B c states are shown in the bottom portion of Fig. H. For any 
given value of /x, the dissociation energy for the B c state is larger than the 
value for the corresponding state in charmonium, and the average radius is 
smaller. 
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Figure 2: Top: B c dissociation energies, bottom: B c radii as functions of inverse screening 
length. 



4 Outlook and Conclusions 

The computation of the ratio Rb c , Eq. (|l|), in the deconfined QGP environ- 
ment poses a challenge, and our estimates yield at present a few percent. In 
any case a very significant, by a factor 100-1000, and observable enhancement 
over elementary processes is to be expected. Based on the known parameters 
of the RHIC we have worked out some interesting numbers shown in table 0. 
The reason that the chemical non-equilibrium for charm is leading to a greater 
yield is the fact that the directly produced charm is here dominating the ther- 
mal production. There is practically no re-annihilation of charmed quarks as 
plasma cools down and hence the significant increase of the tri-lepton channel 
decay yield of the B c . 



Tabic 1: RHIC yields for several heavy quark systems. 



observable 


first year 


design luminosity 


cc-pairs 




2.810 s 


6.5 10 y 


6fo-pairs 




1.2 10 6 


3.2 10 7 


J/* -> 




1.6 10 5 


3.9 10 6 


T(l S )^/x+ M - 




140 


3800 


5 C -> J/* + Zi/ 




no 


QGP 




-Zv 


0.05-0.18 


1.5-4.9 


B c -> J/* + Z^ 




therm+chem. equil. 


QGP 


-+ fj, + fj,~ -i 


- Zi^ 


18 


490 


B c -» J/fy + Iv 




only prim, cc 


QGP(T = 400MeV) 


-> -f 


-Zr/ 


235 


6400 


£ c -> J/* + Zi/ 




only prim, cc 


QGP(T = 200 MeV) 


-> ft + fi~ -f 


- Zi^ 


1000 


27,000 



Our conclusions are: 

(1) The larger binding energy of the B c system and the availability of the 
primary-interaction charmed quarks facilitates production and stability of the 
B c meson; 

(2) The B c should be observable at RHIC (as well as at LHC- Alice), and the 
study of the elementary properties of the B c mesons appears possible; 

(3) The dynamics of the B c system in plasma (given its smaller size and greater 
binding) is different from that of J/tp, and thus B c should produce information 
that is complementary to that produced by J/tp, and T(ls); 

(4) There is potential for a 'smoking gun' evidence of deconfinement indicated 
by a significant enhancement of Rb c , and thus it is important to evaluate this 
quantity theoretically and measure it experimentally in QGP environment. 
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